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Abstract 

Generalizing matrix models, tensor models generate dynamical triangulations in any dimension and support 
a 1/A expansion. Using the intermediate field representation we explicitly rewrite a quartic tensor model as 
a field theory for a fluctuation held around a vacuum state corresponding to the resummation of the entire 
leading order in 1/A (a resummation of the melonic family). We then prove that the critical regime in which 
the continuum limit in the sense of dynamical triangulations is reached is precisely a phase transition in the held 
theory sense for the huctuation held. 


1 Introduction 


Matrix models [T] generate dynamical triangulations and dynamical tessellations in two dimensions. Tensor models 
m II12 0 [6] generalize matrix models and generate dynamical triangulations 018] in any dimension (recently a 
matrix model generating dynamical triangulations is dimension three has been proposed E)- The continuum limit 
of dynamical triangulations is reached when tunning the coupling constants to some critical values [M HU mis]. 
In this regime the number of simplices in the triangulation diverges. Sending simultaneously the volume of the 
individual simplices to zero one obtains a continuous phase of infinitely refined random spaces mm- 

Tensor and matrix models can be regarded as field theories with no kinetic term. The fields are a pair of 
complex conjugated tensors , T^i (where = 1,... A) transforming under the external tensor product 

of D fundamental representations of the unitary group U{N): 
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The action of tensor models is a function of T and T which is invariant under unitary transformations. Like matrix 
models, tensor models are known to possess a 1/A expansion [Ml [151 [161 HZl HI]- The leading order in the 1/A 
expansion for matrices is given by the planar triangulations. For tensor models, the leading order in 1/A is given 
by a class of triangulations called melonic [12] . The melonic triangulations represent topological spheres in any 
dimensions, but one should emphasize that only a subclass of triangulations of the sphere are melonic. In particular 
the family of melonic triangulations is exponentially bounded m- It is not yet known whether the family of general 
triangulations of the sphere in dimension higher than three is exponentially bounded or not [19] . 

In this paper we consider a tensor model with a quartic interaction. We first rewrite this model, using the 
intermediate field representation [201 HI], as a coupled multi matrix model. Examining the vacua of this multi 
matrix model we find for small coupling a unique vacuum state invariant under unitary transformations which we 
call the melonic vacuum. We derive the effective theory for the fluctuations around the melonic vacuum. The 
translation to the melonic vacuum is not a phase transition as no symmetry gets broken (the melonic vacuum 
possesses the full invariance of the theory). We show that this translation completely resums the leading order 
in 1/A: the effective theory for the perturbation fields around the melonic vacuum contributes only to lower 
orders in 1/A. We show furthermore that, at each order in 1/A, only a finite number of graphs of the effective 
theory contribute. Finally, we show that the criticality of the leading order contribution in 1/A, which is the 
dynamical triangulation continuum limit, corresponds precisely to a zero eigenvalue in the effective mass matrix for 
the fluctuation field. 
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We therefore prove that the continuum limit in dynamical triangulations is a phase transition in tensor models. 
This phase transition presumably corresponds to the spontaneous breaking of the unitary invariance of the theory. 
In order to establish this rigorously, one needs in the future to identify the appropriate non unitary invariant vacuum 
state and study perturbations around it. 

The situation is quite different for matrix models. A procedure similar to the one we use here for tensors can 
not be used for matrices: so far at least, the action of no matrix model could be translated to an explicit “planar 
vacuum state” such that the effective theory around it contributes only a finite number of graphs at each order in 
IjN. 

A very interesting question is whether similar results hold for group field theories EDEniinn! and tensor field 
theories [SHESlEelEzlES]. Such theories have a genuine renormalization group flow [291130] and it would be very 
interesting to test whether the phase transition we identify in this paper is associated to an infrared fixed point 
of this flow in which case tensor field theories and group field theories would naturally generate bound states of 
extended geometry. 

This paper is organized as follows. As our results are somewhat technical, we start by discussing in section [2] 
their interpretation. In section [3] we introduce the quartic tensor model and in section |4| we discuss the effective 
theory around an invariant field configuration. In section[5]we introduce the class of Feynman graphs (combinatorial 
maps) relevant for the Feynman expansion of the effective theory which we detail in section [6l Finally in section [3 
we consider the translation to the invariant vacuum state and detail the IjN expansion of the fluctuation field. 


2 Discussion 

The partition functions of tensor and matrix models are generating functions for dynamical triangulations [D Eli- 
Denoting A a (strongly) connected D dimensional triangulation and Top(A) the number of D simplices in A, the 
logarithm of the partition function (the free energy) writes as: 

VF= J_lnZ= ^ 5Top(A)^-(^<^(A) ^ 

-/V A 


where N denotes the size of the matrices (or tensors) and w(A) is a non negative integer. In the case of matrices 
D = 2, a;(A) is the genus of the triangulations [T], while in higher dimensions [HI [1^] a;(A) is the degree. While 
the degree is not a topological invariant, it serves to organize the series defining IF as a series in 1/A. 

In the large N limit only a subclass of triangulations contribute: the planar triangulations in dimension two 
and the melonic triangulations in higher dimension. Both families are exponentially bounded. Similar results hold 
order by order in 1/N and it can be shown that [I] [311133] summing the families of triangulations at fixed order 
in 1/A, the free energy becomes (up to non universal analytic pieces in g which can be disposed of by a suitable 
number of derivatives): 

W -^doigc-gy° + E ■^dqigc-gY'^ ■ 


9>1 


While the critical value gc depends on the details of the model, all the series at fixed order in 1/A become critical 
at the same critical value g^ The (non integer) critical exponents are, to a large extent, universal for a given 

dimension. 

Considering equilateral simplices of volume Vd, the physical volume of a triangulation is proportional to the 
number of D simplices, and the average physical volume is: 


Vol = Vd (Top(A)) = Vd gdglnW ■ 


g-^gc 


Vd 

9-9c 


The continuum limit of dynamical triangulation is reached by tuning g ^ gc and sending at the same time Vd 0 
while keeping the average physical volume fixed. In this regime one obtains a phase of random infinitely refined 
geometries. One can study more involved matrix and tensor models, adding for instance coupling constants for lower 
dimensional simplices. The various continuum phases of the model will then be spanned by these new couplings. 

Dynamical triangulations can also be studied 0133] with no reference to matrix and tensor models. In this 
case one restricts the family of triangulations by some prescription. For instance one can chose to restrict to 
triangulations of the sphere, or to foliated triangulations [33] • Of course one should check that the restricted family 
of triangulations is exponentially bounded (in order to be able to obtain a continuum limit) but, due to technical 
difficulties, this is never done in practice. Still, one can learn numerous lessons from numerical studies. 


2 




We show in this paper that a certain tensor model can be recast as a coupled multi matrix model for D matrices 
with c = 1,... D. We then show that, in this new formulation, the vacuum state of the model is not = 0, 
but it is given at small coupling constant g by the invariant configuration; 


H‘^ = aol 


ao = 


1 - yi - ADg^ 

2Dg 


We translate to the vacuum state and write the partition function of the model as some explicit factor multiplying 
the partition function of the effective theory of the perturbation fields around the vacuum state: 


^(4) = f [dM^] , 
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SiM) = - al) ^ Tr, - -N 
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2 ETrc[M=] -Q{M). 


( 1 ) 


We prove on the one hand that the integral over the perturbation fields is subleading in IfN, namely: 

-W = + ln(l - gDa,) + O (1) , 

and on the other that the mass matrix of the effective theory for the perturbation \m-o eigenvalues 

- Oq) with degeneracy ND'^ - 1 and - Da^) with degeneracy 1. 

Dynamical triangulations continuum limit. The dynamical triangulations continuum limit is obtained from the 
explicit factor in Eq. ©• The leading order free energy writes in terms of g as: 


-Wlo 


^al + \n{l- gDao) 


1 ^ l-^l-4Jg2 

2 ADg"^ \ 2Dg'^ j ’ 


and becomes critical for the critical constant: 


In the critical regime p- = 4^Dg = 1- e one obtains: 


-Wlo 




i+ln2+ie--e3/2 + 0(ey , 

2 2 3 


hence the critical exponent z/q = §• 

Field theory phase transition. The effective field theory for the perturbation field (which is subleading in 
1/A^) undergoes a phase transition when the mass matrix has a zero eigenvalue, that is for I-Doq = 0. At this point 
the invariant vacuum is no longer stable and the symmetry of the theory is broken. This occurs for the coupling 
constant gp,t. solution of: 


1 1 ‘^^9p.t. 

\/~D ‘^hdgpt 


yJl-ADgl,, = {l-ADgl,) 


9p.t. 


1 

2^/D ■ 


We conclude that, in tensor models, the dynamical triangulations continuum limit is a phase transition associated 
to a breaking of the symmetry (unitary symmetry and color permutation) of the model. 

We stress that the phase transition we identify in this paper is of an entirely different nature form the phase 
transitions studied so far in the context of dynamical triangulations and matrix models: 

• Dynamical triangulations. In dynamical triangulations one takes the continuum limit and obtains 

various continuum phases. In Euclidean dynamical triangulations 0 one typically obtains two continuum 
phases (crumpled and branched polymer), while in causal dynamical triangulations [33) one obtains three 
continuum phases. One is then interested in the transition between the continuum phases. These transitions 
have nothing to do with the phase transition from discrete to continuum we identify in this paper. 
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• Matrix models. In matrix models [84] one studies the transition from a one cut resolvent to a multi cut 
resolvent. The phase transition we identify in this paper is very different, and occurs in a completely different 
range of parameters. In matrix models the continuum limit is a change of behavior of the spectral density of 
a one cut solution (one passes from a square root law to a 3/2 power law for the distribution of eigenvalues). 
It is this change of behavior which we show in this paper is a genuine phase transition in the case of tensor 
models. 


3 The quartic melonic model 


A detailed introduction to the general framework of tensor models can be found in m or [35) , and an introduction 
to the quartic melonic model in EQ] and m- An important notion for tensor models is the notion of H-colored 
graph [5]. Such graphs are dual to {D - l)-dimensional triangulations [5]. 

Definition 1. A bipartite edge D-colored graph (D-colored graph for short) B is a graph with vertex set V{B) and 
edge set £{B) such that: 

• the vertex setV(B) is the disjoint union of the set V^(B) of white vertices and the setV^(B) of black vertices, 
with\V'^{B)\ = \V'’{B)\ = k{B). 

• the edge set £{B) is the disjoint union of D sets £^{B) of edges of color c, with c e {1,... D} such that any 
edge connects a black and a white vertex. 


• all the vertices have degree D and all the edges incident to a vertex have distinct colors. 


The fields in a tensor model are a pair of complex conjugated rank D tensors T^ji T„i with = 1... N. 
We call the position c of an index its color and denote the set of colors V = { \ ... D}. We denote a® = (a'^,c e V). 
The tensors T and T transform under the external tensor product of D fundamental representations of the unitary 
group U{N): 


N 

Kn= E 

hi,...6D = l 




N 

E 

6i,...bD = l 




Starting from T and T one can build a class of invariant polynomials called trace invariants. They are in 
bijection with bipartite edge D-colored graphs. The trace invariant associated to B is obtained by associating to 
every white vertex i of B a tensor to every black vertex j of B a complex conjugated tensor Tgp and to every 

_ J 

edge e'^ = {i,j) of color c of 13, a contraction of the indices of color c of the two tensors associated to its end vertices 
i and j: 

(k{B) \ (k{B) 

TrB(T,f) = I n 1 ( n 

The melonic quartic tensor model the probability measure: 


n 


\e'==(i,j)ie(B) 






Z(4)(A) 


2*7r 




n'^'nPmPmP 


( 2 ) 


where, for any set C c T>, we denote i5^c„c = ElceCThe model is stable for g^ < 0. Each invariant T) 

can be represented as a D-colored graph with two white and two black vertices and signifies a gluing of four D- 
simplices, two positively oriented and two negatively oriented. The coupling g counts the number of positively 
oriented simplices. When evaluating the partition function in perturbation theory one obtains Feynman graphs 
having {D + 1) colors. Each such graph is dual to a D-dimensional triangulation. 

The generating function of the moments of is: 

The partition function is the generating function at zero external sources and is denoted by = Z^^\0,0). The 
cumulants (connected moments) of p^^'^ can be computed as the derivatives of the logarithm of Z\ 


'^2k{^aY I TgD , . . . T^-n , Tgl) ) = 


Qi^k) 


dJ^T>dJ^T> . . . dJfjpdJ^T> 


lnZ(4)(J, . 
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The cumulants are linear combinations of trace invariant operators [35]: 


K2fc [ToD , TgC . . . T„I> , TgC ] 


B. k(B)=k e<^ = {i,j)i£{B) " 


where the sum runs over all the D-colored graphs [5] B with k{B) white vertices. Somewhat abusively, we call the 
coefficient R{B,p,^‘^'^) also a cumulant and we define the rescaled cumulant corresponding to the graph B'. 

K(B _ 

) - j^u_2k{B)(D-l)-C(B) ’ 

where C{B) denotes the number of connected components of B. 


3.1 The intermediate field representation 

The intermediate field representation of the quartic melonic models has already been discussed in detail in [20] and 
m- Let us recall here the main steps in its derivation. Each quartic interaction in expressed in the intermediate 
field (Hubbard Stratonovich) representation as an integral over an auxiliary N x N hermitian matrix: 




D-l £ 


3 — y 
2 




f [dH,] 


lAfD-i Tr[H'=H'=]+gN^ 




.,'Ds{o}^Dn{o}I 


,I3n{c}„I3s{c} I 


where the integral over is normalized to 1 for 5 = 0. As there are D quartic melonic interactions, one obtains D 
intermediate matrix fields Each is an N x N matrix and the indices of H'^ have the color c. The integral 
over T and T is now Gaussian and can be computed explicitly to obtain: 






S{H) = y Trdil'=il"]-Tri,[lni?(il)] , 

^ C=1 


where Ti’c is the trace over the index of color c, is the trace over all the indices, and the resolvent is: 

liiH) = --TT > 

where 1 '^ is the identity matrix with indices of color c and for any subset C c D we denote I'' = ( 8 >c 6 C 1 '^- 

In the rest of this paper g will be considered real, positive, and smaller than Observe that in this 

range of the coupling constant the resolvent can be singular. However one can still make sense of the integral 
defining 2'0)( J) by modifying the contour of integration for the field H‘^ in the complex plane so as to avoid the 
singularities. 


3.2 Equations of motion 

Taking into account that: 

d 

———Trx) [lni?(it)] = 5 ^ R{H)g-Dp-D [SpT,-^{c}gT>-.{c}SpcacSqchc) , 

the classical equations of motion in the intermediate held representation write: 

® " or/c ^ -9 y R{H)qDpT> {SpT>^{c}qT>^{c}5pc{,c6qCac) . 

pr>qD 

The = 0 configuration is not a vacuum of the theory as, for 5 0, H‘^ = 0 is not a solution of the classical 

equations of motion. In order to hnd solutions to the classical equations of motion we express them in matrix form: 


H^-g 






= 0 


(3) 
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where Trx)N{c} denotes the trace over the indices with colors in 2? \ {c}. A simple inspection of these equations 
reveals that for g small enough the stable vacuum of the theory is invariant under conjugation by the unitary group 
and invariant under color permutation, = al. The invariant solutions of eq. are then obtained for a satisfying 
the self consistency equation: 

9 


1 - gDa 


This self consistency equation has two solutions: 


The melonic vacuum. A first solution, denoted ao and called the melonic vacuum, is the sum of a power series 
in g: 


ao = 


1 - VI - iDg‘ _ „ 

^o{n+iy.\n / 


We will see below that this solution is the stable vacuum of the theory for g small enough. In the region of 
interest {g eR+, ADg'^ < 1), oq is real, positive, and bounded by: 


ao = 


l-Vl-4Jg2 

2Dg 




1 


1 + ^/U^AD^ Vd 


as 


29 


1 iP ^ increasing with g and attains its maximum for g = 


The instanton. The second solution, denoted ainst, is an instanton solution: 

l + y/l- 4.Dg^ 


^inst ~ 


2gD 


As ooflinst = j), it follows that in the region of interest (g e ]R+, 4,Dg^ < 1), ainst is real, positive, and Oinst ^ 

For g = —4= the two solutions collapse: an| i = i = — 

y 2'/D ^ 2y/D V-^ 

4 Translating the intermediate field 

In order to study the effective theory around a solution of the classical equations of motion we must translate the 
field to = al + This translation is a translation (in the complex plane) by a of the contour of integration 
of all the diagonal entries of H'^. The effective theory is obtained by Taylor expanding in M°. We first discuss a 
translation of H'^ by some arbitrary constant a. We denote M = Y,^=i (gi and we define: 


h{a) = 


9 


I - gDa 


Observe that if g is real and a is real, then b{a) is also real. The melonic vacuum and the instanton are the two 
solutions of the self consistency equation a = b{a). We will use the shorthand notation b instead of b{a) but the 
reader should keep in mind that 6 is a function of a. 

In terms of the perturbation field M, the resolvent and the trace of its logarithm write as: 


i?(al^ + M) = 


(1 - glla)l^ - gAf \l-gDa 


\l-qDa) 1^- 


Ti'i) [lni?(al^ + Af)] = -N^ ln(l - gDa) + Trp 


In 


bM ’ 
I 


I'^-bM 


After translation by a the action becomes: 


I 


D 


1 D 

^ -rD-1 


-N^Da^ + N^-'^a Y, Trc[Ar] + -N^-^ Y Trc[Af‘=Af‘=] + ln(I - gDa) - Trp 


C=1 


C=1 


In 


1 


1'^ -bM 


and the generating function writes in terms of the perturbation field M : 
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D \ „D-1 

YlidM^^Ue ^ 

C=1 / 


Tr4M‘=M‘=]-a Tr^ [M]+Tr^ [in ^ j( ) ,7 


The logarithm of the resolvent contains quadratic terms in M, which must be reabsorbed in the quadratic part 
of the action. Let us define the subtracted vertex function: 


Q{M) = Tit, 


In ■ 


1 


Tiv[aM]-^TiT,^{bMY'^ , 


which, by definition, has no quadratic terms in M. The generating function is then: 


kc=l / 


Integral form of the subtracted vertex function. One can give a closed integral formula for the subtracted 
vertex function Q(M). A Taylor expansion with integral rest yields: 


Trx) 


In 


1 


T^-bM 


= {dt Tit, 


In- 


1 


1® - tbM 
and substituting the derivatives: 


+ -\dtTiT, 

t=0 ^ 


In 


1 


1® - tbM 


t=o 


1 f dt{l-tfdfTiT, 

2 Jo 


In 


1 


1® - tbM 


a: Tit, 


In 


1 




= (r-l)!Trc 


{I'^-tbMy 


:{bMy 


we obtain: 


r ^ 

Q(M) = (b-a)TiT,[M]+ dt(l-tyTiT, 

Jo 




{bMf 


that is Q{M) is the sum of a linear piece and a piece starting at order M^. 


Linear term. At the price of modifying the covariance, the translation of H'^ by a generic a introduces some 
univalent vertices in the theory. Choosing a = ag comes to choosing the weight of the univalent vertex as the 
counterterm of the melonic graphs HHS], and by translating to the melonic vacuum, one subtracts the contribution 
of the entire melonic family: as ag = b{ag), the linear term (&-a)Trx)[M] which represents the univalent vertices 
vanishes. As expected, in this case = 0 is a solution of the classical equations of motion due to the fact that the 
remaining piece of Q{M) starts at order M^. 

No such interpretation exists for the translation to the instanton Oinst although the linear term cancels also in 
that case. 

Free energy. We can compute the large N free energy of the model in the perturbative regime using Schwinger 
Dyson equations. According to m, the large N covariance K 2 and the large N free energy are: 


99 1 - ag 

l-K2 + Dg^Kl = Q ^ K2= - 2 =" 


= ^im lnZ(4)(0,0)) = 1 + ln(iL2) -K 2 + ^KI = -^KI + hiiK^) . 
The large N free energy writes in term of ag as: 


Woo = “^“0 “ ln(l - gDag) 


that is precisely the explicit prefactor in Eq. (|31) for a = ag. We will re-derive this result below. 
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It ensues that the translation to the melonic vacuum H‘^ = (IqI resums the contribution of the melonic sector, 
and the remaining integral over M (which, as we will see below, only contributes at lower orders in 1/N) is precisely 
the effective theory of a perturbation field around the melonic phase. 

We note that the case of matrices, D = 2, is quite different. A careful inspection shows that for matrices, even 
after translating to a = oq, the integral over the perturbation M contributes at leading order in IjN. This is 
natural because for matrices the melonic graphs are only a subset of the planar graphs, and one cannot subtract 
the contribution of the planar sector by a single counterterm. 


4.1 The effective theory for the fluctuation field 

The space of Hermitian NxN matrices is a N'^ dimensional real vector space Ti with inner product {A\B) = Tr[Ai?]. 
The quadratic part of the action is a quadratic form on the direct sum of D such Hilbert spaces 0^i Denoting 
a generic element in this direct sum space by: 

( M^\ 


M = 




the inner product in 0^^ is: 


D 


{P\M) = . 

C=1 


The effective theory for the fluctuation field M is somewhat involved due to the fact that, while the quadratic 
part of the action is a quadratic form on the direct sum space the interaction Q{M) is a trace in the 

tensor product space (8>^i This is the root of all the subtleties we will need to deal with below. 

We now work out the effective covariance for the fluctuation field. The quadratic part of the action of the 
fluctuation field in Eq. o is: 

D 

j^D-i ^ XrdM'=M'=] - TrviibMf] = 

C=1 

D D 2 / ^ 

= iV'°"i(l-6^)^Trd(M‘=)2] + 7V'^~262^(TrdM‘=]) - ^ Tr^M^] . 

c=l c=l '' \c=l / 

We denote by Ic the element of 0^^ consisting in the identity matrix in and zero on the other components 
and by / the element consisting in the identity matrix in all the 'H'^s: 


Ir = 


/o\ 

1 

VO/ 


I = 


^l\ D 

i -ZP- 

a/ 


The orthogonal projectors on the one dimensional subspaces generated by Ic and I are: 


P„ M = - (MIP. M) - 


N 


PM = 


{ic\ic) 

Tr,[M-] ~ I ^ Tiy[A/T 


N 


mi) 


■1= 


DN 


'-i^{M\PM) = 


DN 


||;^Tr,,[Ar']j . 


Let us denote I the identity operator on 0^i The quadratic part writes then as (M | O M) where the operator 
O is: 

O = -b^)l + N^-^b^ E Pc - N°-^Db^ P . 

0=1 

The covariance of the effective theory is the operator In order to compute it, it is most convenient to 

diagonalize O. 












Lemma 1. The operator O is diagonalized as: 


0 = 


D 


i-EPc 


+ N 


D-1 


D 


EPc-P 


+ N^-\l-Db^)V , 


where I - Pc; Pc “ P and P are the projeetors on the eigenspaces of O. The eigenvalues of O and their 
degeneracies are: 


Ai = 


, dim |lm I - ^ Pc j 


= DN^ - D 


A 2 = ^ , dim [ Im 


D 


EPc-p 


\ Lc=i 

As = N°-\l - Db'^) , dim(Im[P]) = 1 . 


= D-1 


Proof. It is enough to check that the operators I - Pc, Z£:i Pc - P and P are three (mutually orthogonal) 
orthogonal projections. The three operators are obviously Hermitian and: 

Pc;PcM = , 


N 


N 


(|:pc)p^ 

p(ePc)m=E 




E' 


DN ctl N 


-L = PM 


^ Trc[M=] (i|/c) 


N DN 

The rest of the proof is a straightforward computation: 


J = PM. 


D 

i-EPc 


D 

i-EPc 


D 


EPc-p 


D 


EPc-p 


D D D 

= I-2^Pc+ E PcPc-I-EPc; 

c=l c,c'=l c=l 

Z PcPc'-p(|:Pc)-(ePc)p + p 

,c'=l \c=l / \c=l / 


c,c- 
D 


= EPc-p, 


p2 = P 


D 

i-EPc 


D 1 ^ ^ / ^ \ 

^Pc-P =EPc- E PcPc'-p+ EPc p = o 

c=l . c=l c,c'=l \c=l / 


D 


i-EPc 


c=l 


P = P 


D 


EPc-p 


p = 


(I/-) 

(I 


p = o 


p^ = o. 


Concerning the dimensions of the images of each of these operators, by definition P is one dimensional, Pc and 
Pc' are orthogonal c+ c' hence [l - Pc] is DN^ - D dimensional and finally the sum of the three projectors 
is the identity. 

□ 


It is now easy to compute the covariance of the effective theory by inverting O: 


= 


I 


ND-i(i-b2) 
1 


D 

I-EPc 


C=1 


1 


ND-1 


D 

EPc-p 


.c=l 


1 


-I- 


7V-d-i(1-62) 


(I'-l 


ND-i{l-Db'^) 
Db^ 


ND-i(l-Db‘^) 


P . 
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Observing that the normalized Gaussian integral with covariance C for any field cj) can be written as a differential 
operator |5S] : 


/ 


[#] 

TditC 




1^=0 ’ 


and, taking into account the normalization of the Gaussian integral, the generating function in Eq. o becomes: 

-Ar^(-2fH.+ln(l-g_Da)) 






1° %*) (l-BCa) 


M‘^=0 


(5) 


Using the explicit form of the projectors Pc and P, the differential operator in the above equation writes in terms 
of derivatives with respect to M‘^ as: 


/ 


\dM 


O 1 




1 


D 


dMl c=i 

62 


ETU 


_9_^ 

dM^ dM^ 

D 

E Trc 


D 


7V^(l-62)^Vi IdM' 


ETU 


d 


Trc 


d 


dM^ 




d 


dM^ 


Trc 


d 


dM<^' 


The effective theory can be signihcantly simplihed. The crucial observation is that the differential operator 
( 9 ^ I functions which depend only on M = Y,c=i yve observe that the action of an 

operator: 


Trc 


r 1 / ^ \” n- 1 +...n 

L9M^J\c=i / 


nl 
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® ... 0 = n y iVr 0 




n-l 


D\no 


does not depend on the color c. It follows that all the terms Trc[^^] can be replaced by the symmetrized operator 
^ Trc[g^] hence the operator | can in fact be replaced with: 


(A 

\dM dMl 


D 




ETrc 


d d 

dM^ dM^ 






]f-ETrc[^l] . 

)\Dti ldM-\) 


The new operator, O ^, which defines the covariance of the effective theory can be written in terms of I and P 
only: 

1 ^ ^ 1 1 


0“" = 


-1 + 


7 V£>-i(i_ 52 ) iV^-i(i- 62 )(l-£) 62 ) lV-D-i(l- 62 ) 

Gorrespondingly, the quadratic part of the action becomes: 


(I-P) + 


7V^-i(l-D62)' 


D f) _ \ ! D \ 2 

= 7V'^“1(1 - 6^) ^ Trc -6^ ^ ■ 

c=l T) \c=l / 

With this remark, the generating function can then be rewritten in the simpler form: 

7V°(-efi+l„(l-gDa)) 


V'(l-62)i5Af"-^(l-£)62) 

(1-1,2) ivfi(\-”2°(iEl,2) ('P 2i£.i 2:°iTrc[^^]) 


[ 1_L 

g 2 Ar-D-1( 


^Qi^)+ ffD-1 a-gDa)j{ 


M'==0 
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or, in integral form: 




/ 




(7) 


where the integral over M° is normalized to 1 for g = 0,b = 0, J = J = 0. 

We conclude that, after a translation by a of the intermediate field, the quartic tensor model becomes (up to an 
overall factor) a field theory for the fluctuation matrix fields with effective action: 


S(M) = -N 


D-1 


D 


(1-6^) ^Trc 


iiV^-2£z. 
2 D 


Tr,[M=] -Q{M). 


Stability 

The stability of the effective theory is encoded in the effective mass matrix of the fluctuation field, which is: 


dMdM Im=o 


= N^-^{l-b^){l-V) + N^-'^{l-Db^)V . 


The effective mass matrix has then two eigenvalues: ^{1 - b^) with degeneracy ND^ - 1 and ^(1 - Db^) 

with degeneracy 1. Let us analyze this effective mass matrix for the melonic vacuum and the instanton. 


The melonic vacuum. The melonic vacuum oq is a solution of the equation 6(ao) = uo- Furthermore, oq 
is bounded from above by for g < and attains this value exactly a.t g = It follows that all the 

eigenvalues of the effective mass matrix around the melonic vacuum are positive for g < and at precisely 

g = one of these eigenvalues becomes 0. The melonic vacuum is therefore the stable vacuum of the theory for 
g < and at g = the effective mass matrix around this vacuum develops a zero eigenvalue, hence the effective 
theory for the fluctuation field undergoes a phase transition. This phase transition corresponds to a breaking of the 
full symmetry of the model (conjugation by D unitary transformations and color permutation). 

The instanton. The instanton is also a solution of the equation b(ao) = oq, but this time oq > The effective 
mass matrix around the instanton has at least a negative eigenvalue in the region g < , hence in this region the 

instanton is unstable. 


5 Edge multicolored maps 

In this section we introduce the notion of edge multicolored maps. We will show in the next section that the 
Feynman expansion for the fluctuation field M is indexed by these objects. 

Definition 2. A combinatorial map with external cilia is: 

• a finite set of half edges S which is the disjoint union of the set of internal half edges iSint and the set of 
external cilia iSext- We denote the elements of S by he S. 

• a permutation a on S. 

• an involution a on the set of internal half edges iSint having no fixed points. The involution a on iSint is extended 
to an involution on S (which we denote also by a) with hxed points by imposing that a{h) = h, V/i e 5ext) 
that is by imposing that all the external cilia are fixed points of a. 

The cycles of the permutation a are the vertices of the map: the half edge cr{h) is the successor of the half edge 
h when turning around a vertex of the map. These maps have a well defined notion of faces. 

Definition 3. The cycles of the permutation aa fall in two categories: 
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• the cycles of aa which do not contain any ciliuin, that is such that: 

^ ‘^ext 5 

for any h in the cycle are called the internal faces of the map. 

• the cycles of aa which contain cilia. They are further subdivided into sequences of half edges separating two 
consecutive cilia: 


h, aa{h),..., (aa)’^(h) 


such that • 


e tSext 

VO < s < r, a{aaY{h) 
a{cray{h) € 5ext 


i Se-xt 


called external strands. 

Definition 4. The corners of a combinatorial map are the couples: 

(/i,cr(/i)), yheS . 


The corners can be pictured as the pieces of vertices comprised between two consecutive half edges or between 
a half edge and a ciliuin. 

An edge multicolored map with cilia is a map with cilia whose edges are furthermore colored by a subset C of 
colors, C c {!,... D} (in the next section we will encounter only maps such that C consists in a unique color, but 
for now we keep the discussion general). 


Definition 5. An edge multicolored combinatorial map with external cilia A4 is: 


a finite set S which is the disjoint union of the sets 
of external cilia.' 


of internal half edges of the colors C and the set iSext 



• a permutation a on S. 

• for every C cT), an involution on the set of internal half edges of colors C, , having no fixed points. 

We extend the involutions to the whole of S by setting a^^Yh) = h, VheS \ . 

An example of an edge multicolored map with external cilia is presented in Figure [TJ 

c 



Figure 1: An edge multicolored map with external cilia (the cilia are the half edges in the figure). 


Let us consider an edge multicolored map Ai. For any color c e {1,... D}, by erasing the edges of colors C with 
c ^ C, one obtains the submap Ai‘^ with external cilia and scars of A4. 

Definition 6. The submap with external cilia and scars of the edge multicolored map with cilia: 

M = ((|J 4?)U‘5ext,cT,{a(‘')|c C 2?}j , 

is (see Figure\^for an example): 
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• the finite set S‘^ of half edges of color c which is the disjoint union of: 

— the internal half edges of color c: 

sLt= U 4? 

C, C€C 

— the external cilia; 

5ext = 5ext 

— the scars; 

cc =11 

^scars I_I ^int 

C. etc 

that is: 

5“ = 5r„,U5s^carsU5exf 

• the permutation a on S. 


• the involution on 5“^^ having no fixed points defined by: 

a^{h) = (h) if h e and ceC . 

The involution oT is extended to an involution on S with fixed points by imposing that each scar and each 
cilium is a fixed point of a‘^, i.e. a'^{h) = h for he S‘^ \ iSj'Jjf 


c 



Cl 



vCl 


Cl 


C2 * c, Cl 


Figure 2: Subniaps of color c, Ci and respectively C 2 c, Ci of the map in Figure [T] where we did not represent the 
scars in the various submaps. 


The map can be disconnected, and some of its connected components can consist in an isolated vertex. The 
corners of an edge multicolored map A4 and its submaps are still defined as the couples (h,cr(/i)), for all the 
half edges h. 

Definition 7. The faces of an edge multicolored map with cilia and scars are the cycles of the permutation 
aa'^. They are divided into 

• internal faces. They are the cycles of aoT which do not contain any external cilium, that is, for every h in 
the cycle: 

The internal faces can however contain any number of scars, 


• external faces. They are the cycles of tra^ containing cilia. As before, they are further subdivided in seguences 
of half edges separating two consecutive cilia: 


h,aa^{h),, (tTa°)’’(/i) such that 


■ Vs, 0 < s < r, a°(CTa^)®(/i) i 
a^(aaT(/i)e5=xt 


called external strands. 
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If the map A4‘^ has a connected component consisting in an isolated vertex, it counts as a face. 

The internal faces and the external strands of color c of an edge multicolored map M are the internal faces and 
the external strands of the submap We denote by the number of internal faces of color c of M and 

by ^lnt(-Ad) the total number of internal faces of A4. 

Remark that the faces of an edge multicolored map are colored by a unigue color. For example, the faces of 
colors c, Cl and C2 + c,ci of the map in Figure [T] are the faces of its submaps represented in Figure [2l Also, observe 
that the scars play no role in the definition of the faces. Their only relevance is to keep track of the corners of the 
map, and in particular to ensure that an edge multicolored map M and its submaps have the same corners. 

Definition 8. The boundary graph dAi of an edge multicolored map with cilia A4 is the edge D-colored graph 
obtained by associating a black and a white external vertex to each of its cilia, and connecting the white vertex 
associated to the cilium h with the black vertex associated to the cilium h' by and edge of color c if there exists an 
external strand of color c of the map Ai (i.e. of its submap A4‘^) going from h to h'. 

An edge multi colored map can be represented either, as we have done so far, as a map with colored edges or as 
a multi ribbon graph. This representation is obtained as follows (see Figure [3] for the representation of the map in 



Figure 3: Multi ribbon graph representation of an edge multicolored map with external cilia. 


Figured]); 

• we represent every internal vertex by D concentric strands colored 1 to D. By convention we orient the 
strands clockwise. 

• we represent every ciliated vertex as D strands connecting a white and a black external vertex. By convention 
we orient the strands from the white to the black vertex (one can further collapse the white and the black 
vertex into a cilium). 

• for any edge of colors C we connect the strands with colors in the set C on its two end vertices by a ribbon. 
The edges of the graph are thus made of \C\ parallel ribbons. 

The advantage of this second representation is that the faces and the external strands are easily read off: the 
internal faces are the closed strands, while the external strands are the strands connecting a white and a black 
external vertex. The boundary graph of Ai is the graph obtained by erasing the internal faces of Ai in this 
representation. 

It is important to note that the maps we obtain in the next section never have more than one cilium per ciliated 
vertex. 


6 Feynman rules 

We will now discuss the Feynman rules for the effective theory defined by the Eq. (|6|) and 0. 
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6.1 The effective vertices 

The generating function J) can be evaluated perturbatively. Using for instance the differential operator 

form for the Gaussian integral, Eq. (jb)), and expanding the exponential we obtain: 






X Z 

tlint ittext^O 


1 

^int !^ext! 


g 2 \ QM 


dM I 


r ntlint / 

[QiM)] ^ 


7V-C-1 


1 

(l-gDa) 



As usual such an expression is evaluated in terms of Feynman graphs. The graphs are obtained as follows. We 
represent the terms Q{M) as internal vertices and the JJ terms as external vertices. On the external vertices we 
add a cilium signaling the presence of the external sources Jj. The derivatives in the operator | O 

on the vertices, and each derivative creates a half edge. The operator connects pairs of half edges into edges. 
We now evaluate the action of the derivative operators on the various vertices. For the ciliated vertices we have: 


dM 


a<=b'= 


1 

-b y 

1 


1 


[I'^-bMi 

q^p-n n^mn \ 


q^nn 


mPpn j 


( 6 . 




) • 


For the internal, non ciliated vertices, we need to estimate the derivatives of the subtracted vertex function: 


Q{M) = Trp 


In ■ 


■Trp[aM]-iTr7,[(6M)2] . 


Clearly, starting with three derivatives, the last two terms do not contribute. For the lower orders we use: 
d 


dM 


2 


Trp [aM] = (a T^M^ ® = a ^ ^ 

V c=l / pT:>qD 

D D \ 

^Tri,[(M^)2® ^ Tri,[Ar ®M'='® l®x{c}x{c'}] 


^ J 2 dM^ 


t)C \ C=1 


c,c' = l 
c^c' 


= ('^P®x{c},jl3s{c}l5pc:a'=^g'=b‘=) , 

p-Dg-D 


dM 


■Tr 


V 




In 


1 


I'^-bM 


-b E 


1 


. 1® — bM iq'DpTt 


((5p-Ds{e}q-Ds{e}(5pC£,el5qCbc) . 


It follows that the first derivative of the vertex function is: 


d 


dM^ 


-[q{m)]= Y, 


-bM 


■ai^ -b^M 


(5pr>N{c},jr>N{c}(5pC(jC(5qCf,c) , 


and the second derivative is computed using: 


d 


dM 


a’^b'^ 

= b^ 


I'^-bM 


-al^ -b^M 




E ( 

'PmP \ 



1 


1 

l 


q'Dri'D 



[l ]m^p^ I (^n'^''{c>^'Dx{c} ) . 


6.2 Broken and unbroken edges 

When evaluating the Gaussian integral in Eq. dH), the two terms in the differential operator | corre¬ 

spond to two kinds of edges: 

• unbroken edges colored by a color c coming from the term: 


Tr, 


\ d a ] 
.aM^ dM<=_ 
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• broken edges coming from the term: 



The broken edges have no color. 

We draw the broken edges as dashed edges and the unbroken edges as solid edges having a color. As usual, as 
the fields are matrices, the order of the half edges around a vertex is relevant and the Feynman graphs are in fact 
combinatorial maps (a.k.a. ribbon graphs [I]). An example of a map with broken and unbroken edges is presented 
in Figure m 


i, 

/ s 
/ ^ 



Figure 4: Map with cilia, broken and unbroken colored edges. 


In the perturbative Feynman expansion, the generating function writes as a sum over maps A4 with ciliated 
vertices whose edges fall in two categories: the unbroken edges f„(A4), which have a color c, and the broken edges, 
Sbi-M.)- Let us denote the successor permutation of the map A4, and let us denote the corner by 

[/i]. Estimating the derivatives acting on the vertices, we obtain for every corner of the map an associated matrix 
element of an operator: 

• for any corner [h] of a ciliated vertex or of a non ciliated vertex of coordination at least three we obtain: 

1 






for a corner [h] of a non ciliated vertex of coordination one we obtain: 

1 


1^-bM b 


[fe] [fc] 


for the couple of corners [h] and belonging to a non ciliated vertex of coordination two we obtain: 



1 


1 

1 

[I'^-bMi 

CL^ b'^ 

[hVih] 

[l^-bM] 


- [1 ]^T) J-D [1 ] -D J-D 

,-D k-d -'“[M [M 




The edges contract the indices of the operators incident at their ends in a specific pattern, according to the 
nature of the edge (unbroken or broken). Extracting the logarithm we obtain a sum over connected maps. We 
denote nint(A4), next(A4), Eu{M) and Eb{M) the numbers of internal vertices, external vertices, unbroken edges 
and broken edges of M. Recall that S{M) denotes the set of all the halfedges of A4, and iSext(A4) the set of cilia 
of A4. With thees notations, the perturbative expansion of the generating function of the cumulants writes as: 


+ ln(l - gDa) 


)- 


ln(l-6^) -^\n{l-Db^)+ 


Z 


^2E^{M)+2Eb(M) 


E 


M connected ’^int (A4 ) Iriext (A4 ) ! 


\hiS(M) ' ’ ^ V (1 gDa) ^ 
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(Sb‘: 


x[ n b^D{D-l) (, . 

where we have used the notation: 

• for any corner [h] of a ciliated vertex or of a non ciliated vertex of coordination at least three: 

1 




Jm=o 




l^-bM 


[M [M 


for a corner [h] of a non ciliated vertex of coordination one: 

1 




1® - bM b 


-^1^-bM 


[M [»>] 


for the couple of corners [h] and \<JM{bi)\ belonging to a non ciliated vertex of coordination two: 

i®1 _ _ It'D 


[<*] [fe] [»■('*)] 


1 


1 

[I'^-bMi 

n'^ h'^ 

[l^-bM] 


- [1 ]jjI3 f,© [1 jjj-D (,13 

H hT> -■“[(>] [(>] -'VCh)] 


7 Translating to the melonic vacuum 

If we translate to the melonic vacuum, several significant simplifications can be made. Using &(ao) = oq and taking 
into account that M is set to zero we obtain: 

• for any corner [h] of a ciliated vertex or of a non ciliated vertex of coordination at least three: 


^(0)„ 




= 1 


V 

aJi.b'F, ’ 


[M W [M [M 

for a corner [h] of a non ciliated vertex of coordination one: 

b{ao)-ao-,v n 

-b(U)a® b® = -77-7-1„I3 u-D = U , 

^“wdM b{ao) 

for the couple of corners [h] and [o'/ 4 (/i)] belonging to a non ciliated vertex of coordination two: 


3^(0)a® b® ^(0)o® b® 

[fe] [h] [o-Cf*)] l<r(.h)] 


]a?,b^,[^ hp,,,,bP,,,, ]a^,b^,[^ hp,.,,bP,,,, 

[fc] [fc] [o-(fi)] [h] [fc] [o-Cfi)] [cfCli)] 


Thus only maps whose internal vertices have coordination at least three (which we denote Ai^) survive. More¬ 
over, a map Ai^ has unbroken colored edges and broken edges. Erasing all the broken edges leads to a map whose 
edges are colored by a color c = {!,... D}. We denote this map by Ai^ (as it is formed only by the unbroken edges). 
Remark that Ai^ can be disconnected, even though Ai^ is connected (see Figure [5]). 

As Ail^ is a usual edge colored map it has a well defined notion of boundary graph dA4^ and (internal and 
external) faces. The indices are identified along the faces of As the numbers of unbroken edges of Ai^ and 
Ail^ are equal, and the number of ciliated vertices of A4 is the number of white vertices of we obtain: 


lnZ('‘)( J, J)=-N^ + ln(l - gDa) j - ^ ln(l - 6^) - ^ ln(l - Dh'^) + 

y Trains (J,J) 1 

connected k{dA4iy. nint(Al^)! 


2S^ 

Oq 


1 


r a^DiD-l) 1 

Eb(M^) 

(1- 


. 1 - “o. 
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Figure 5: The maps and 


where k{dAi^) denotes the number of white vertices of the boundary graph (see definition [51) and 

k(dMl) 

Tra^3(J, J)= n ^afJaf n ■ 

i=l e‘^=(vi,Vj )€£{Ai^) 

The rescaled cumulants are computed by taking derivatives and dividing by the appropriate power of N: 




1 


2Eu(M^)+4:Eb(M^) 


[!)(£» 


connected 

8A13=S 


nint(M3)! (1 - 

^ jY~DEb(M^)-(D-l)E^(M^)-D+(D-l)k(dMl)+C(dMl)+Fi„t(Ml) 


( 8 ) 


7.1 The 1/N expansion for the fluctuation fleld 

In this section we show that the cumulants admit a 1/N expansion. We label the connected components of 

and we denote C(Nl^) their number. 

Definition 9. For each connected component we define its deficit ®'S' 

= 2 + E^iMl^p)) + {D - 2)V{Ml.^p))- 

-{D- mOMl^^p^) - CidMl^p)) - FUMl,^p)) . 

Furthermore, we define the excess of broken edges of , Lb(At^), and the excess of unbroken edges o/ 
as: 

L,{M^) E E,{M^) - C{Ml) + 1 , L,,iMl^^p^) e E^(Ml^^p^) - V(Ml.^p^) + 1 . 

Theorem 1 (1/A expansion). The rescaled cumulants of expressed as the formal series: 

1 p2E^{M^)riEb{M^) [D(D - 

ni„t(Al3)! (l-gDao)'=(^^-) {1 - {I - 

dM^.B 

^ j^-DLbiM^)-(D~2) . (9) 

The right hand side of Eq. (|9]) is a series in 1/N. More importantly, there are only a finite number of maps At^ 
contributing to any fixed order in 1/N. 

Proof. A straightforward computation leads from Eq. ([5]) to Eq. (O taking into account that the boundary graph 
and the internal faces of Al^ are distributed among the connected components and that the numbers of 

vertices, ciliated vertices and unbroken edges of Af^ and Al^ are equal. 

As Ai^ is connected, all the connected components of Al^ must be connected in between them by broken edges, 
hence Lb(Ai^) > 0. Similarly, as is connected, Lu{M./^.^p^^) > 0. 

In appendix [A] we prove the following lemma: 


A/t3 connected 
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Lemma 2. The number of internal faces of a connected edge colored map with k(dM) cilia is bounded by: 


FintiM) <1-{D- l)k{dM) - C(dM) + (D- l)V(M) + [E(M) - V(M) + 1] . 

As a consequence of this lemma, for any connected component its deficit is non negative, 

^(■^u-(p)) - li' follows that the right hand side of Eq. ([9]) is a series in IjN. 

In order to prove the second part of the theorem, we note that, as D > D - 2, the scaling in N is bounded by: 

r C(Ml) 

p=i 

As the ciliated vertices in the map have coordination at least one, while the non ciliated vertices have 
coordination at least three, we have: 

2[Eb{M^) + Eu{M^)] > 3 [V{M^) - k{dMl)] + HdMl) ^ V{M^) < \[Eb{M^) + E^iM^) + k{dMl)] , 

hence the scaling with A^ of a term on the right hand side of Eq. (|S]) is bounded by: 

(E + 2k{dMl) 

^ 3 ' 

The number of external vertices of is fixed, as it equals the number of the white vertices of the invariant B 
whose cumulant we evaluate. It follows that at any order vnljN only connected maps with at most a finite number 
of edges contribute, and the second part of the theorem follows. 

□ 


Setting the external sources to zero we obtain the free energy of the quartic melonic model: 


InZ(^) =-Af-° 
+ N^ 


Da^ 


+ ln(l - gDa) 


ln(l ln(l - Db^)+ 


connected ^itltCAd^)! 

3Ad3.0 


2 E^(M^)+ 2 Ei,(M^) 

X Oq 

1 

1 - a, 



r 


^oD(D-l) 




The sum over Ad^ does not contribute at leading oder in N. Indeed, a term in the sum with Lb(Ad^) = = 

0 would be a tree, hence would have univalent vertices, which is impossible as all the vertices Ad^ are at least three 
valent. 


A Proof of lemma [2] 

In this appendix we show that the number of internal faces of a connected edge colored map Ad cilia is bounded 
by: 


Fi„t(Ad) <!-(£»- l)fc(9Ad) - C{dM) + {D- l)E(Ad) + [A(Ad) - E(Ad) + 1] . 

Proof. In order to establish this bound, we first introduce a new colored map Ad* associated to Ad which will allow 
us to keep track of the boundary DM while modifying Ad. Starting from Ad, for each cilium h, we introduce 2D 
half edges on the corner bearing the cilium h, one at its left and one at its right for each color c € 2? (that is, for 
each color c, one half edge precedes h and the other succeeds h when turning around the vertex). We then connect 
these new half edges into dashed edges following the edges of i9Ad: if the external strand of color c starting at the 
cilium h ends at the cilium h', we connect the half edge of color c following h with the half edge of color c preceding 
h'. This construction is represented in figure [3 
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Figure 6: A map A4 with two cilia. For simplicity, only the colors 1 and 2 have been represented. 


Each of the external strands of M corresponds to an edge in the boundary graph dM, hence A4 has exactly 
Dk{dAi) external strands. By construction M.* has a new edge for each external strand of M. which closes the 
external strand of M into an internal face of It follows that: 

FintiM*) = Fint(M) + Dk{dM) . 

On the other hand, as the new dashed edges exactly follow the external strands of AI, it ensues that dAi* = dA4. 



Figure 7: The map M* obtained from the map M in Fig. |6]by adding the dashed edges. 


The reason to introduce A^* is the following. Below we will delete solid edges of A^* (that is edges which 
belonged to AI) and track the change in the number of internal faces of AI* under these deletions. The crucial 
point is that these deletions modify the structure of the internal faces of Ad* but do not modify its external strands, 
hence the boundary graph dM* will remain unchanged. The introduction of Ad* allows us to modify a map while 
keeping its boundary graph unchanged. 

Let T be a tree in Ad and let us denote T* the graph made of T and all dashed edges of Ad*. We have: 

dT* = dM* = dM . 

The map T* can be obtained from Ad* by deleting E{M) - V(M) + 1 solid edges. As there is only one face 
going trough any of the solid edges, the deletion of such an edge in Ad* can either increase of decrease the number 
of faces by one, hence: 

FintiM) + DkidM) = FintiM*) < FintiT*) + EiM) - ViM) + 1 . 

We call a leaf of T* a vertex which becomes univalent when one erases all the dashed edges. A moments 
reflexion reveals that the leafs of T* are the univalent vertices of T. In order to complete the proof of lemma O it 
is now enough to show that for any T* we have: 

FintiT*) + CidT*) < 1 + kidT*) + iD- l)Vir*) . 

This bound is obtained by tracking the evolution of Fint(T*) +C(9T*) under the iterative deletion of the leaves. 
This deletion is somewhat involved, as it can change the boundary graph. 

We say that an external strand of T* is looped at the eilium h if it starts and ends at the cilium h. On the 
boundary graph dT*, this corresponds to an edge connecting the black and white vertices associated to h. In hgure 
[7] we can identify two looped external strands. Let us denote I a leaf of T* and let us denote T* the tree obtained 
from T* after deleting i as follows. 

Deleting a non ciliated leaf. If t has no cilium, deleting it simply means deleting the vertex I and the edge 
connecting it to the rest of T*. If the edge connecting I to the rest of T* has color c, the deletion of I erases all 
the faces of colors T> \ {c} running trough 1. The boundary graph is unaltered by this deletion, therefore: 

FintiT*) = FintiT*) + (D-1), Cidrr*) = Ci&t*) ^ 

FintiT) + CidT*) = FintiT) + Cidff*) + (D - 1) . 
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Deleting a ciliated leaf. If I has a cilium h, deleting it consists in : 

• Step 1. For all the external strands starting or ending at h that are not looped, we cut the corresponding 
dashed edges into half edges. We then reconnect the dashed half edges into edges the other way around, 
respecting the colors. This is represented in figure [8] below. Observe that after performing this step, all the 
external strands going trough i are looped. 



Figure 8: The first step of the deletion of a ciliated leaf : the dashed edges have been reconnected in such way that 
all the external strands are looped at the cilium h. 


• Step 2. The leaf now has only looped external strands and it is connected to the rest of the T* by a solid 
edge only. The cilium h represents a connected component of the boundary graph consisting in a black and 
a white vertex connected by D edges. We erase the vertex £, its cilium, h, and the edge connecting i to the 
rest of the graph, as in figure |9] 



Figure 9: The second step of the deletion of a ciliated leaf. 


Boundary graph. Upon deleting a ciliated leaf, the boundary graph changes: dT* + dT* ■ We have two cases, 
each with two sub cases: 

• None of the external strands of T* are looped at h. One needs to apply Step 1 for all the colors and then 
Step 2. There are two sub cases: 

— the black and white vertices associated to h in dT* belong to the same connected component of dT* ■ 
Then, in the boundary graph. Step 1 creates at least a new connected component, and Step 2 deletes 
exactly one connected component. Thus: 


C{dT*)<C{df*) . 

— the black and white vertices associated to h in dT* belong to different connected components of dT* ■ 
Then, in the boundary graph. Step 1 can not decrease the number of connected components, and Step 
2 deletes exactly one connected component. Thus: 

c{dT*)<c{df*) + i. 

• At least one external strands of T* is looped at h. There are two sub cases: 
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— not all the external strands are looped at h. The black and white vertices associated to h in dT* belong 
to the same connected component of dT* ■ One must apply Step 1 at least once, and then step 2. As 
before, Step 1 creates at least a new connected component, and Step 2 deletes exactly one connected 
component, hence: 

C{dr*)<C{df*) . 

— all the external strands are looped at h. The black and white vertices associated to h in dT* belong 
to two different connected components of i9T*, and one must apply Step 2 directly. This decrease the 
number of connect components of the boundary graph by 1: 

C{dT*) = C{df*) + l . 

Internal faces. Let us denote e'^ the solid edge of color c connecting i to the rest of T* ■ We have several cases: 

• the external strand of color c is looped at h. Then there is only one internal face of color c in T* running 
through e°, which can not be erased by deleting £, hence: 

F-ntiT*) = Ff^,{f*) . 


• the external strand of color c is not looped at h. Then there are either one or two internal faces of color c in 
T* running trough e°, hence the number of internal faces of color c can not decrease by more than one: 

PUT*) < FUf*) + 1 . 


• the external strand of color c' + c is looped. Then there is only one internal face of color c' through i, which 
is erased: 

fL{T*) = Fi,{f*) + 1 . 

• the external strand of color d + cis not looped. Then there is just one internal face of color c' through £, 
which is not erased: 

fL{T*) = Fi,{J*) . 

Combining the counting of the connected components of the boundary graphs with the counting of the internal 
faces we obtain four cases: 

• no external strand is looped at h. Then: 

Aint(r*) + G{dT*) < [Fi„t(r*) +1] + \c{dft*) +1] < Fi„t(r*) + c{df*) + d . 


• all the external strands are looped at h. Then: 

^i„t(r*) + G{dT*) < [Li„t(r*) +D-1] + [G{df*) + 1] < Fintif*) + G{df*) + D . 

• the external strand of color c is looped at h, but at least one external strand of color c' c is not looped at 
h. Then: 

Pint {T*)+G{dT*)< [Fint {f*)+D-2] + [G{df*)]<Fint{f*) + G{df*)+D . 

• the external strand of color c is not looped at h, but some external strands of colors c't c are looped at h. 
Then 

FintiT*) + G{dT*) < [Tntif*) +D]+ [G{df*)] < Fintif*) + G{df*) + D . 

In all cases, by deleting a ciliated leaf: 

Ff^,{T*) + G{dT*) < F,l,{f*) + G{df*) + D . 

Iterating up to the last vertex, we either end up with a vertex with no cilium or with a ciliated vertex with D 
looped external strands. Counting the number of internal faces and connected component of the two possible final 
graphs we conclude. 

□ 
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